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ABSTRACT 

We have computed stellar evolutionary models for stars in a mass range 
characteristic of Cepheid variables (3 < m/M Q < 12) for different metallicities 
representative of the Galaxy and the Magellanic Clouds populations. The 
stellar evolution calculations are coupled to a linear non adiabatic stability 
analysis to get self-consistent mass-period-luminosity relations. The period - 
luminosity relation as a function of metallicity is analysed and compared to the 
recent EROS observations in the Magellanic Clouds. The models reproduce the 
observed width of the instability strips for the SMC and LMC. We determine a 
statistical P-L relationship, taking into account the evolutionary timescales and 
a mass distribution given by a Salpeter mass function. Excellent agreement is 
found with the SMC PL relationship determined by Sasselov et al. (1997). The 
models reproduce the change of slope in the P-L relationship near P ~ 2.5 days 
discovered recently by the EROS collaboration (Bauer 1997; Bauer et al. 1998) 
and thus explain this feature in term of stellar evolution. Some discrepancy, 
however, remains for the LMC Cepheids. 

The models are also in good agreement with Beat Cepheids observed by the 
MACHO and EROS collaborations. We show that most of the 1H/2H Beat 
Cepheids have not yet ignited central helium burning; they are just evolving off 
the Main Sequence toward the red giant branch. 
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1. Introduction 

Over the past five years, tremendous efforts have been devoted to the search for 
dark matter through microlensing effects. This systematic search has provided two large 
databases of variable stars in the Magellanic Clouds due respectively to the EROS group 
(Renault et al. 1996, Beaulieu & Sasselov 1996 and references therein) and the MACHO 
group (Welch et al. 1996 and references therein). The wealth of data collected for Cepheids 
by both projects provides a unique tool to analyse period - luminosity relations and their 
sensitivity to metal abundances. A significant number of Beat Cepheids, oscillating in both 
the fundamental (F) and first overtone (1H) modes or the first and second (2H) overtone 
modes have also been observed in the Clouds. 

The analysis of the EROS Cepheids (Sasselov et al. 1997, S97) in the small (SMC) 
and Large (LMC) Magellanic Clouds has recently revived the old debate regarding the 
metallicity dependence of the period - luminosity (PL) relationship (see Tanvir 1997 for a 
review). A careful comparison of the observed PL relations in both clouds enabled S97 to 
disentangle from differential distance and reddening a significant metallicity effect, which is 
not predicted by theoretical analysis (Chiosi, Wood & Capitano 1993, CWC93). 

Analysis of Beat Cepheids based on the OPAL opacity data set (Iglesias & Rogers 
1991; Iglesias and Rogers 1996) yield now good agreement for galactic Bump and Beat 
Cepheids (Moskalik et al. 1992; Christensen-Dalsgaard 1993, Christensen-Dalsgaard 
and Petersen 1995), but discrepancies still remain regarding the Magellanic Clouds 
(Christensen-Dalsgaard & Petersen 1995; Buchler et al. 1996; Morgan & Welch 1997, 
MW97). 

None of the afore-mentioned analysis is really consistent since the stability analysis are 
based on envelope models using input physics different from the ones entering the stellar 
evolution calculations. Moreover, the width of the instability strip and in particular the 
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position of the blue edge are deduced from stability analysis, varying arbitrarily T cS for a 
given mass - luminosity (ML) relationship, independently of any stellar evolution constraint 
on T e ff. The aim of this letter is to present fully consistent stellar and pulsation calculations 
based on the most recent opacity data set. In §2 we present the stellar and pulsation 
calculations. Section 3 is devoted to the comparison with the observed PL relationships for 
fundamental mode pulsators in the SMC and LMC. The Beat Cepheids are analysed in 
section 4 and conclusions are given in §5. 

2. Stellar evolution and pulsation calculations 

Our stellar evolution calculations include the updated OPAL opacities (Iglesias & 
Rogers 1996) for T > Q000K. For lower temperatures, we use the Alexander & Fergusson 
(1994) opacities. Convection is based on the mixing- length theory with a mixing length 
parameter l mix = 1.5 H p . The onset of convective instability is based on the Schwarzschild 
criterion, without either overshooting or semi-convection prescriptions. 

The linear non-adiabatic stability analysis is based on a radial pulsation code originally 
developed by Lee (1985; Lee and Saio 1989). The convective flux F conv is assumed to be 
frozen in, i.e. the perturbations 5F conv are neglected into the linearized energy equation. 

The stellar evolution and pulsation codes are coupled. This yields a fully consistent 
mass - luminosity - period - evolutionary time relationship. We have computed stellar 
models in the mass range 3-12 M & with various chemical composition (Z, Y) = (0.02, 
0.28), (0.01, 0.25) and (0.004, 0.25), representative of respectively the galactic, LMC and 
SMC environments. The detailed evolutionary tracks will be published elsewhere and we 
focus in the present letter on the main results. Figure 1 shows the evolutionary tracks of 
Z=0.004 stars in a theoretical Hertzsprung-Russell diagram (HRD), for m = 3 - 7 M Q . 
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Unstable models are indicated by either open circles or crosses. The distinction is purely 
artificial: the crosses refer to the usual Cepheid "blue loop" instability which occurs during 
the central helium burning phase while the open circles correspond to the first rapid crossing 
of the HRD when the stars evolve off the Main Sequence toward the red giant branch on 
a thermal timescale. This phase, of typically ~ 10 4 years, is referred to in the following as 
the "first crossing instability". Since this phase is about 1/100 shorter than the blue loop 
one, the large data base now available may allow the detection of some Cepheids in the first 
crossing instability. 

As shown in Fig. 1, the extension of the 3 M blue loop (crosses, solid line) is 
significantly reduced, compared to higher mass stars. In this case, the blue edge is 
determined by the turn- over of the evolutionary track. This reduction of the blue loop 
extension is expected from stellar evolution since blue loops during core helium burning 
start at a minimum mass m min and vanish for m >1OM , depending on the initial 
composition and the convective mixing treatment (cf. El Eid 1994). Note that m min 
decreases with metallicity. We find m m i n ~ 3M for Z=0.004, m min ~ 3.75M for Z=0.01 
and m min ~ 4.5M m for Z=0.02. 



3. Period - luminosity relationships 

In Fig. 2a-b we compare our calculations with observed fundamental mode pulsators 
in the LMC and SMC in a P - My diagram. The EROS (B E , R E ) magnitudes were 
transformed into the Johnson- V magnitude according to Beaulieu et al. (1995). We adopt 
the extinction E(B-V) = 0.10 for the LMC and E(B-V) = 0.125 for the SMC with R v = 
3.3 and distance moduli (m — M) = 18.5 for the LMC and (m — M) = 19.13 for the 
SMC (cf. Beaulieu et al. 1997a). Comparison is also made with the Laney & Stobie (1994, 
LS94) data, which extend to longer periods than the EROS data. The transformation of 
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the theoretical quantities (L, T eff ) into My are based on the Allard and Hauschildt (1998) 
most recent atmosphere models. These models originally developed for M-dwarfs do not 
extend below log g < 3.5. We therefore use the bolometric corrections at constant log g = 
3.5 but take metallicity effects into account. We verified however that gravity effects in the 
V-band are small. 

We determine a theoretical statistical P - My relationship under the form 
My = A logP + B with a weighted least square fit by minimizing the quantity : 

Q = J2^(M Vl - A\ogP t -Bf 

i 

where the summation extends over all the stellar masses. The coefficients a^m^t) depend 
on the mass distribution and on the evolutionary times. We adopt a Salpeter mass function 
(MF) dN/ dm oc m~ 2 35 and the time dependence derives directly from the coupled evolution 
and pulsation calculations. Although the present day MF in the MC may differ from an 
initial Salpeter MF, we verified that variations of the slope of the MF between -2 and -4 
barely affect the slope of the PL relationship. Such a steep MF favors the lowest mass stars 
because (i) of the number of stars itself and (ii) of the longer time spent in the instability 
strip as the mass decreases. The slope of the PL relationship is thus hardly affected by 
stars with m > 7M . 

In Fig. 2a, comparison is made between the SMC observations and the Z=0.004 
models with masses m= 3-12 M Q . The first crossing unstable models are indicated by 
open circles but only for m < 4M Q . Although included in the calculations, the first crossing 
instability phase for m > 5M is statistically insignificant, since the time spent during this 
phase is more than 300 times smaller than the time spent by a 4 M G in the instability 
strip (first crossing and blue loop). We thus predict the existence of Cepheids with periods 
P < 1 day, as the signature of the first instability strip, but the observation of these objects 
requires larger statistics. The models agree reasonably well with the observed width of the 
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instability strip, although they do not reach the observed blue edge of the EROS data. 
The faintest objects (M v > —1.5, logP < 0.4) seem to indicate a change of slope in the 
PL relationship, becoming steeper for log P < 0.4, as discovered in the EROS-2 Cepheid 
sample and analysed carefully by Bauer (1997) and Bauer et al. (1998). We note that such 
a trend is observed in the theoretical relation near ~3 M . This change of slope is thus 
real and stems from stellar evolution, illustrating the reduction of the He blue loop as mass 
decreases (cf. §2). Note also that the period of the minimum mass undergoing a blue loop 
(~ 3M ) is consistent with the faintest observed SMC Cepheids. Finally the average slope 
of our P - My relationship is A = -2.92, in excellent agreement with the slope derived by 
S97. 

Fig. 2b shows the results for the LMC with Z=0.01 and masses from 3.75 to 12 M . 
As in Fig. 2a, the models reproduce correctly the observed width of the instability strip and 
the overall agreement with observations is good. We note that the minimum theoretical 
"unstable" mass undergoing a blue loop m min ~ 3.75M for Z=0.01 does not correspond to 
the faintest objects observed, which correspond to ~ 4.25M . We predict that fundamental 
pulsators with logP<0.3 days should be in the first crossing instability. The slope of the 
PL relationship with a minimum mass of 4.25M is A = -2.50, shallower than the one 
observed in the LMC by S97. In order to test the influence of chemical composition, we 
have recomputed the whole grid of models with (Z=0.008, Y=0.25) and (Z=0.01, Y=0.28), 
without any substantial modification of the theoretical slope. For solar metallicity models, 
we find also a slope A = -2.55 shallower than the one observed in the Milky Way. A detailed 
analysis of the possible sources for such discrepancies is under progress. 
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4. Beat Cepheids 

Figure 3a-b display the Petersen diagrams for SMC (Fig. 3a), LMC and galactic (Fig. 
3b) Beat Cepheids. The data are taken from EROS (Beaulieu et al. 1997b) and MACHO 
(Alcock et al. 1995; Alcock et al. 1997) collaborations. Comparison is also made with 
galactic F/1H pulsators (Szabados, 1988). Agreement with the F/1H mode pulsators is 
excellent for the three metallicities, except for the shortest periods of the LMC objects. The 
first crossing unstable models are also displayed, but predict slightly too high 1H/F period 
ratios and extremely short fundamental periods. A variation of Y from 0.25 to 0.28 does 
not affect significantly the period ratios and the Z=0.008 models give slightly higher period 
ratios than observed in the LMC, by 1-2%. The present models however yield a better 
agreement with observations than the one obtained by MW97 using the Chiosi (1990) ML 
relationship with overshooting. 

The most striking and new result shown in Fig. 3 concerns the 1H/2H pulsators, which 
are reproduced by models in the first crossing instability. Indeed, low luminosities are 
required to reproduce the observed period ratios at very short periods, as noted by MW97. 
The luminosities during the He burning blue loop are already too large and yield too high 
Pi, as indicated by the crosses in Fig. 3, as already noted by MW97. In order to reproduce 
the observed shortest periods, they need masses m < 2 M Q . Fig. 1 shows however that the 
evolution during the first crossing takes place at a significantly fainter luminosity and allows 
to reach the shortest periods observed with masses m > 3M Q . We find that at least the 
four faintest 1H/2H Beat Cepheids of the Beaulieu et al. (1997b) sample can be explained 
by models in the first crossing instability. This is statistically reasonable, considering that 
the total number of single mode Cepheids (F or 1H) found by EROS in the SMC is 400. 
Since the lowest masses contribute dominantly to this total number, the probability to have 
4 Beat Cepheids in this sample in the first crossing instability i.e 4/400 ~ 1% is consistent 
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with the ratio of the time spent in the first crossing instability phase to the total time spent 
in the instability strip (F and 1H instability strips included) (cf. §2). 

5. Conclusions 

We have computed self-consistent calculations between stellar evolution and linear 
stability analysis for SMC, LMC and galactic chemical compositions. This first consistent 
analysis yields good general agreement with the observed P - My diagrams in the SMC and 
LMC and reproduces reasonably well the width of the instability strip in both cases. 

The theoretical slope of the PL relationship is in excellent agreement with the empirical 
one derived by S97 for the SMC. The puzzling trend of a change of slope observed for log 
P < 0.4 and discovered by Bauer (1997) is reproduced by the models and thus can be 
explained by purely evolutionary effects due to the reduced extension of the blue loop with 
decreasing mass. 

For the LMC, the theoretical slope is ~ 15% less steep than the observed one. 
The minimum mass found to undergo a blue loop (M = 3.75M , Puueioop ~ 2 days) 
corresponds to periods shorter by 20% than the ones of the faintest observed LMC Cepheids 
(Pobs ~ 2.5 days). Variations of the metallicity and initial helium abundance do not 
solve this discrepancy. This shows the limitation of the present calculations and points 
out a remaining puzzle for the LMC Cepheids. Investigation of remaining uncertainties 
in the current stellar evolution models is under progress but requires numerous lengthy 
calculations. We indeed note that overshooting can change the minimum mass undergoing 
a blue loop and thus the minimum period. However, the good agreement found for the 
SMC seems to indicate that the present standard evolutionary models are already basically 
correct. 
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Regarding Beat Cepheids, we show that 1H/2H pulsators with the shortest periods 
can only be reproduced by stars in the first crossing instability phase, with luminosity 
significantly lower than the one in the usual blue loop phase. This new result is supported 
by statistical arguments (cf. §4). 

Finally, the present letter shows that consistent calculations between stellar evolution 
and pulsation analysis may give the right answer to some unexplained features of Cepheids. 

We are grateful to the EROS group and in particular to F. Bauer for discussing the 
EROS-2 Cepheid sample prior to publication. We are also grateful to U. Lee for providing 
the original pulsation code and to W. Glatzel for helpful discussions. The calculations have 
been performed on the cray T3E of the Centre d'Etudes Nucleaires de Grenoble. 
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Figure Captions 

Fig. 1- Theoretical H-R diagram for stars with Z=0.004, Y=0.25 and different masses. Open 
circles and crosses indicate the location of fundamental unstable modes. The distinction 
between both symbols is explained in the text (§2). 

Fig. 2a.- Period- magnitude relationship for the SMC. Observations are from EROS (plus) 
and Laney & Stobie (1994, crosses). The models are for Z=0.004 and Y=0.25. Dots 
correspond to the core He burning unstable models. Open circles correspond to unstable 
models in the first crossing instability phase (see text). The straight line is the average 
statistical P - My relation derived from the present models (see text). The slope and 
minimum mass are indicated in the upper left corner. 

Fig. 2b.- Same as Fig. 2a for LMC data and models with Z=0.01 and Y=0.25 

Fig. 3a.- Petersen diagram for SMC models with Z=0.004. Observations are from MACHO 
(full circles) and EROS (full squares). Crosses correspond to unstable models during the He 
burning blue loop. Open circles correspond to unstable models in the first crossing instability 
(see text). The corresponding masses are indicated for both F/1H (lower sequence) and 
1H/2H (upper sequence) mode pulsators. 

Fig. 3b.- Same as Fig. 3a for LMC models with Z=0.01 and Y=0.25. The F/1H galactic 
pulsators of Szabados (1988) are added (full diamonds) and compared to Z=0.02 models 
(the lowest sequence, +). The corresponding masses are indicated for Z=0.01. The Z=0.02 
models are only shown during the He blue loop and correspond to 4.5 to 6 M from left to 
right. 
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